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' Abstract. We study the weighted heat trace asymptotics of an operator of 

Laplace type with Dirichlet boundary conditions where the weight function 
exhibits radial blowup. We give formulas for the first few terms in the expansion 
(-H 1 in terms of geometrical data. 

3 ' 
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1. Introduction 

1.1. Motivation. The asymptotic analysis of the heat trace provides a natural 
link between the spectrum of Laplace type operators T> acting on functions on a m 
dimensional Riemannian manifold M and the underlying geometry of M. For small 
. time t it links the distribution of the large energy part of the spectrum of T> to 

local geometric invariants of M and its boundary which show up in its asymptotic 
expansion. These invariants play an important in many physical phenomena e.g. 
in quantum statistical mechanics when taking the large volume limits or in the 
Casimir effect [9]. Typically the coefficient of the leading £ -m / 2 term in the heat 
trace expansion for small t is determined by the interior (volume) of M. In many 
situations a detailed study of the boundary behaviour of the heat kernel associated 
with dt + ~D is desirable. One way of obtaining this information is putting a weight 
£ — . in the evaluation of the heat trace. In the setting of the heat content of M this 

corresponds to giving M a non uniform specific heat. It is well known that the 
\q . diagonal element at (a;, x; t) of the Dirichlet heat kernel associated to e~ tT> vanishes 

' like r 2 where r = r(x) is the distance to the boundary. This allows the weights 

to diverge like r~ a , where Re(a) < 3. We will show, using pseudo differential 
calculus, that a modified asymptotic series still exists in this case. For example if 
1 < a < 3 the leading behaviour of the heat trace is fO-™-")/ 2 with a coefficient 
determined by an integral over the boundary of M. 



> 
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^ 1.2. The heat equation. We adopt the Einstein convention and sum over re- 

peated indices. Let M be a compact smooth Riemannian manifold of dimension m 
and with smooth boundary dM. Let V be a smooth vector bundle over M and let D 
be an operator of Laplace type on the space C°°(V) of smooth sections to V. This 
means that the leading symbol of T> is given by the metric tensor or, equivalcntly, 
that we may express in any system of local coordinates X — (^X . ...j X ) and relative 
to any local frame for V the operator T> in the form: 

(La) V = - {gT Idd^ + Ald Xv + A } . 

In Equation (|l.a|) . we let 1 < fx, v < m, we let A\ and Aq be smooth endomorphisms 
(matrices), and we let be the inverse of the metric g^ v := g(d x , d Xv ). Note 
that the Riemannian measure dx on M is given by 



dx := gdx 1 ...dx m where g :— *J det(<7 M „) . 
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Thus, for example, the scalar Laplacian A^ := 5d is of Laplace type since 

(1-b) Am = - {g^d Xv d^ + g^d^ {ggT} &„) . 

We shall use the Dirichlet realization of the operator T>. For t > and for 
4> 6 L 2 (V), the heat equation 

{d t + V)u{x; t) = 0, u(-;*)|eM = 0, limu(-; t) = <£(■) in L 2 (^), 

has a solution u = e~ W (f) which is smooth in (x; t). The operator e~ tT> has a kernel 
Px>(x,x;t) which is smooth in (x,x\t) such that 

u(x;t) = / pv(x, x; t)<j>(x)dx . 
Jm 

In the case of the scalar Laplacian Am, there is a complete orthonormal basis {4>i} 
for L 2 (M) where the € C°°(M) satisfy <j>i\dM = and A M <j>i = h<t>i- The 
corresponding Dirichlet heat kernel pm '■= PA M is given in terms of the spectral 
resolution {c/)j,Ai} via 

p M (x, a;; t) = y] e~ tA * ^ (x)^j (x) . 
/ 

1.3. Heat trace asymptotics in the smooth setting. We use the geodesic flow 
defined by the unit inward normal vector field to define a diffcomorphism for some 
e > between the collar C e := dM x [0, e] and a neighborhood of the boundary in 
M which identifies dM x {0} with dM; the curves r — ► (yo, r) for r € [0, e] are then 
unit speed geodesies perpendicular to the boundary and r is the geodesic distance 
to the boundary. Let F <E C°°(M) be an auxiliary weight function which is used 
for localization. On C e , expand F in a Taylor series 

oo 

F(y, r)~J2 F *(yy Where F i = h( d rY F \ r =0 ' 
i=0 

Henceforth we shall let Tr denote the fiber trace and Tr^2 denote the global L 2 
trace. We then have: 

(l.c) Tr L 2(Fe- tv ) = [ F(x) Tr{p v (x, x; t)}dx . 

Jm 

We note for future reference that on the diagonal, the heat kernel pgm (x, x; t) for 
M. m and the heat kernel ph(x, x; t) on the half space H := {x : x\ > 0} of the scalar 
Laplacian are given by 

(l.d) pR m (x,x;t) = {Airty m/2 and p H (x, x; t) = (4vrt)- m/2 (l - e~ r2/t ) . 

Let dy be the Riemannian measure on the boundary. To simplify future expressions, 
we set 

2{F} = [ Fdx and l bd {F} = [ Fdy . 

JM JdM 

We will also use the notation I{Fdi>}[U] when it is necessary to specify the domain 
of integration U and/or the measure dv. 

Theorem 1.1. Let M be a compact smooth Riemannian manifold. Let T> be the 
Dirichlet realization of an operator of Laplace type. Let F G C°° (M) . 

(1) There is a complete asymptotic expansion as t J. of the form: 

Tr L 2 (Fe- tv ) - t~ m ' 2 J t n a n (F, V) + r^' 1 ^ 2 J t e / 2 a b e d (F, V) . 

n=0 e=o 
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(2) There are local invariants a n = a n (x,T>) defined on M and there are local 
invariants af\ = a b e d i (y,'D) defined on dM for < i < £ so that 

i 

a n (F 7 V)=l{Fa n } and af (F, V) = ^ Z^^oft} . 

i=0 

We refer to [HI H3] for a proof of Theorem 11.11 where more general results are 
obtained in the context of elliptic operator theory and elliptic boundary conditions. 
We shall illustrate these formulas in Theorem 11.31 below. We add a caution that the 
notation we have chosen differs slightly from what is employed elsewhere. 

1.4. The Bochner Laplacian. Before discussing the formulas in Theorem 11.11 in 
further detail, we must introduce the formalism of the Bochner Laplacian which 
will permit us to work in a tensorial and coordinate free fashion. If V is a con- 
nection on V , then we use V and the Levi-Civita connection defined by the metric 
to covariantly differentiate tensors of all types. Let ';' denote the components of 
multiple covariant differentiation - in particular, 4>-^ v are the components of V 2 <^. 
If E is an auxiliary endomorphism of V, we define the associated modified Bochner 
Laplacian by setting 

(l.e) V(g,V,E)(f>:=-g^(f> ;vft -E ( l ) . 

Let L ' )iva and Y ^ be the Christoffel symbols of the Levi-Civita connection on 
M. We have, adopting the notation of Equations (jl.ap and (|l.e|) . the following (see 

my- 

Lemma 1.2. If T> is an operator of Laplace type, then there exists a unique con- 
nection \7onV and a unique endomorphism E on V so that T> = T>(g, V, E). The 
associated connection 1-form u) of V = V(2?) and the associated endomorphism 
E = E{T>) are given by: 

(1) u ll = ^{jg ltu Al+ S r e T„ lt Id). 

(2) E = A - g^id^w^+u^ - lo.T^). 

1.5. Formulas for the heat trace asymptotics in the smooth setting. Let 

indices i, j, k, I range from 1 to m and index a local orthonormal frame {ei, e m } 
for the tangent bundle of M. Let Riju be the components of the Ricmann curvature 
tensor; our sign convention is chosen so that -R1221 = +1 on the sphere of radius 
1 in IR 3 . On the collar C £ , we normalize the choice of the local frame by requiring 
that e m = d r is the inward unit geodesic normal. We let indices a, b, c, d range 
from 1 through m — 1 and index the restricted orthonormal frame {ei, e m _i} for 
the tangent bundle of the boundary. Let L a b '■= g{^ e a Cb,e m ) be the components 
of the second fundamental form. The scalar invariant L aa is the unnormalizcd 
mean curvature (i.e. the geodesic curvature) and will play a central role in our 
investigation. One has the following formulas; note that F2 = \F^ nm : 

Theorem 1.3. Let M be a compact smooth Riemannian manifold. Let T> be the 
Dirichlet realization of an operator of Laplace type. Let F € C°°{M). 

(1) a (F,V) = (47r)-'"/ 2 I{Tr(FId)}. 

(2) ax{F,V) = i(47r)-™/ 2 Z{Tr(6F£; + FR l3Jl Id)}. 

(3) a b d (F,V) = -i(4^)-(" i - 1 )/ 2 X M {Tr(F Id)}. 

(4) a\ d {F,V) = i(4^)- m /2 X 6<i{ Tr ( 2j F 0j L aa id-3F 1 Id)}. 

(5) a bd (F,V) = -3i i (47r)-(™- 1 )/ 2 X M {Tr(F (96£;+ [16^ - 8R amma 

+7L aa L bb - 10L ab L ab ] Id) - 30FiL aa Id +48F 2 Id)}. 

Formulas for the invariants a n (F, T>) and a b e d (F, V) are known for n,£ = 2, 3, 4, 5. 
We refer to [S] for further details as the literature is vast and beyond the scope of 
the present paper to survey. 
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1.6. Singular weight functions. Fix aeC. Let F be a smooth function on the 
interior of M such that Fr a £ C°°(C e ); the parameter a controls the growth (if 
Re(a) > 0) or decay (if Re(a) < 0) of F near the boundary, assuming that Fr a 
docs not vanish identically on the boundary. We may expand F \c e in a modified 
Taylor series: 



oo 

F(y,r) ~J2 Fi (y) ri ~ a whcrc F *W = h( d rY(r a F)\ r=0 . 

i=0 

1.7. Geometry near the boundary. The Riemannian measure is in general not 
product near the boundary, i.e. dx ^ drdy, and this plays an important role in 
our development. Let indices er, g range from 1 to to — 1 and index the coordinate 
frame {d Vl , d Vm _ 1 } for the tangent bundle of the boundary. One may express the 
metric on the collar C £ in the form: 

ds li = 9<yQ(y, r ) d y a ° d y e + dr 2 . 

Fix j/o G dM and choose the local coordinates so that <7<xe(j/o,0) = S ag . Then we 
have that: 

Lag = g{d r ,Va ma d Xe ) = T ag m = -\d r g ag , 

(l.f) g M {y , r) = y/dct {Id +d r g ag (y 0l 0) • r + 0(r 2 )} = 1 - rL aa (y ) + 0(r 2 ), 

dx = (1 - rL aa )drdy + 0{r 2 ) . 

Example 1.4. Let x\ = Ccos6* and X2 = Csm# be the usual polar coordinates on 
the unit disk in M 2 . One then has that ds 2 = d( 2 + ( 2 d6 2 so dx = CdOdC,. The 
geodesic distance to the boundary circle is given by r = 1 — Q; thus gee = (1 — r) 2 
and L aa = 1 so dx = (1 — r)drdy. 

1.8. Regularization. Before discussing the asymptotic expansion of the heat trace 
in the singular case, we must first discuss regularization; an analogous regularization 
was required when discussing the heat content for singular initial temperatures in 
@]. Let H be smooth on the interior of M with Hr a £ C°°(C £ ). Then 

dx = (1 - rL aa )drdy + 0(r 2 ), 

Hdx = {H a r- a + (#! - H L aa )r 1 - a }drdy + 0{r 2 ' a drdy). 

For Re(a) < 3, define: 

lR Cg {ff} :=l{Hdx}[M -C s ] 

+l{Hdx - [H r- a + {Hi - H L aa y- a ] drdy} [C s ] 

(l.g) +l bd {H }x{ Trs; ^ 

y ln(ej ii a = 1. 



- F L aa | X / 2-a lf Q ^ 2 ' 

X ' 1 \ ln(e) if a = 2. 



This is independent of the parameter e and agrees with I{H} if Re (a) < 1. Because 
the integrand over C e is 0(r 2 ~ Rc ^) and Re(2 — a) > —1, 2p,, cg is well defined. 

The regularization T^ e& is a meromorphic function of a with simple poles at 
a = 1,2. At these exceptional values, 2p,, cg is defined as the constant term in 
the appropriate Laurent expansion, thus dropping the pole. We shall apply this 
regularization to functions of the form H(x) = F(x)a n (x,T>). 

1.9. Heat trace asymptotics in the singular setting. The Dirichlet heat kernel 
satisfies pv(x, (y, f), i)|f=o — 0. Since px> is smooth for t > and C e is compact, we 
may use the Taylor series expansion of px> to derive the estimate: 

\pv(x,(y,f);t)\<C(t)f on C e . 
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A similar estimate holds for \pv({y, r),x; t)\. We set x — (y, r) to derive the estimate 
on the diagonal: 

(l.h) \pv((y,r),(v,r);t)\ <C(t)r 2 on C e . 

Thus if Re(a) < 3, then Equation Ql.cp shows that Tr L 2 (Fe - * 13 ) is convergent. 

We shall begin our investigation in Section [2] by establishing the following result 
by a direct computation as this motivates our entire investigation; note that only 
limited smoothness is required of the boundary in this result. It has the further 
advantage of confirming by completely different means some of the constants that 
will be computed again in Sections 4 and 5. Let C here and elsewhere denote Euler's 
constant. 



X(r) 



Theorem 1.5. Let M CM 2 be an open, bounded, and connected planar set with C 
boundary. Let < Sq < e. Set F(x) := Fo(y)r~ a x(r) where \ <= C°°(C e ) satisfies: 

1 if 0<r<e o , 
if s<r. 

(1) 7/0 < a < 1 and t J. 0, then: 

Tr L *(Fe- tA ") = ^- t {T{F} - \T (i=a) ^"^/^{Fo} 
+ 4^S) r (t 2 ) t^ 2 l bd {F L aa }} + 0(1) . 

(2) If a = 1 and t 1 , then: 

Tr i2 (F e - tA -) = ^ {jReg{i^} - | Ht) ■ 2 bd {Fo} + f l bd {F }} + 0(i"*) . 

This result extends to a very general setting. The following Theorem generalizes 
Theorem 11.11 to the singular setting where, in contrast to Theorem 11.51 we assume 
the boundary is C°°. We also refer to [4] for further details where an analogous 
result was proved for the heat content asymptotics. In Section [3J we will use the 
pseudo-differential calculus to show that: 

Theorem 1.6. Let M be a compact smooth Riemannian manifold. Let T> be the 
Dirichlet realization of an operator of Laplace type. Let a n = a n (x,T>) be the interior 
local heat trace asymptotics of Theorem 1.1. Fix a G C with Re (a) < 3. Let F be 
a smooth function on the interior of M such that Fr a € C°°(C e ). 

(1) J/a ^ 1, 2, there is a complete asymptotic expansion as t J, of the form: 



oo 



Tr L 2(Fe- tv ) ~ r m ' 2 J2 t n l Rcg {Fa n } + r^ 1 ^ 2 i ( *~ a)/2 < a (F, V) . 

71=0 1=0 

(2) If a = 1, 2, there is a complete asymptotic expansion as t J. of the form: 

oo oo 

Ty L 2(Fe- tv ) ~ t~ m ' 2 £ t n l Rcg {Fa n } + t _ (" 1-1 )/ 2 ^ t^~ a ^ 2 a\ d a (F, V) 

n=0 £=0 
oo 

+r m ' 2 ln(t) ^ t k / 2 a b k d a {F, V) 



fc=0 



(3) There exist local invariants a\ a i = a\ a i (y,2?) on <9M, which are holomor- 
phic in a for a =/= 1, 2, so i/iai 



i=0 

TTie invariants ah • have simple poles at z = 1,2 and 



e,z,i 

if a = 1,2. 



r ,bd _ I bd r> 00 M 
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(4) The ln(t) coefficients in Assertion (2) are given by: 

( -\l bd {{Fa n \ ) } ifk = 2n and a 

a bd a {F,V) = \ -^{(Fa^i - (Fa n ) L aa } if k = 2n and a 
[ p = 2n + l. 

Throughout this paper, let 

Ka := 2^ (h2^) ' 

The boundary invariants for a ^ 1,2 and for £ = 0, 1, 2 are given by: 

Theorem 1.7. J/a/ 1,2, then one has: 

(1) a bd a (F,V) = KQ (4^)~ m / 2 Z bd {Tr(-F Id)}. 

(2) a\ d a {F,V) = KQ _ 1 (4 7 r)-" 1 / 2 X M {Tr(-F 1 Id + 2^i yJ F 1 L aa Id)}. 

(3) a\ d a {F,V) = Ka _ 2 (A7T)- m / 2 l bd {TT(-F 2 ld+^^F 1 L aa Id 

+^FoR amma Id — 8 ^_ 7 6 ) FoL aa Li,b Id + 4 ^ t _ 5 6 - ) FoL a bL a b Id 

We remark that we recover Theorem 11.31 by setting a = in Theorem 11.71 We 
omit details as the calculation is entirely elementary. The boundary invariants for 
a = 1, 2, and for I = 0, 1, 2 are given by: 

Theorem 1.8. 

(1) Mien a = 1, 

(a) < X (F,£>) = (47r)"™/ 2 2: M {Tr(f F Id)}. 

(b) a^(F,£>) = (4^)- m /2^iM {Tr( _ j p iId+ |^ oiaaId)} . 

(c) a^(F,P) = (470-™^ {Tr (-±F 2 Id +±FL aa Id +±R amma Id 
~-^L aa Li,b ld+j^L a bL a b ld+j^Rijji Id + |E) } . 

(2) Mien a = 2, 

(a) a b d 2 (F,V) = (47r)- m / 2 v^l bd {Tr(F a ld)}. 

(b) a\%(F,V) = (4 7 r)-™/ 2 2: brf {Tr(fF 1 Id-[f + ±]F L aa Id)}. 

(c) < 2 (F,2?) = (4^)- m / 2 ^ M {Tr(-iF 2 Id+|F 1 L aa Id 

+F)( J^Ramma — ^[L aa Lbb + -^L a (,i a fc + \R%jji ) Id +F F)} . 

Here is a brief guide to the remainder of this paper. In Section [2j we will make 
a special case calculation to establish Theorem 11.51 A probabilistic estimate of R. 
Lang [10] and of H. R. Lerche and D. Siegmund [11] plays a central role. In Section 
[3] we shall use dimensional analysis (scaling arguments) and various functorial prop- 
erties to study the heat trace invariants. We will derive Theorem 11.61 (4) from the 
asymptotic series in Theorem [T76] (3); another derivation will be given subsequently 
in Section [5] We shall examine the general form of the invariants and establish the 
following result. 

Lemma 1.9. 

(1) There exist universal constants {R a , k^, k^, k^, k^} so that: 

(a) a b d a (F,V) = (4ir)- m / 2 I bd {Tr{-R a F Id)}. 

(b) a bd a (F,V) = (4^)- m / 2 X M {Tr(-S Q _ 1 F 1 Id+4F L Qa Id)}. 

(c) a bd a (F,V) = (4^)- m /2 X brf{ Tr (_ Sct _ 2F2ld+K i_ ijFlLaaId 

+F [nl,R amma + K^L aa L b b + K 5 a L a bL ab \ ld-K a F Q [E + \Rijji Id])}. 

(2) Ifa^ 1,2, then R a = K a and k\ = (^) 2 (a-3) ■ 

(3) Ri = §. 
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In Section |U wc evaluate the remaining universal coefficients of Lemma ll.9l using 
the calculus of pseudo-differential operators and complete the proof of Theorem 1 1.71 
by showing: 

Lemma 1.10. Adopt the notation of Lemma \1.9[ If a =/= 1,2, then: 

K 3 _ a-l p ( l-a \ .A _ 7-8a+a 2 n I l-ct \ 5 _ 6a-5-a 2 p f l-a 'i 

^ 24 1 V 2 / ' — 32(a-6) 1 V 2 / ' ~ 16(a-6) 1 V 2 J ' 

We conclude the paper in Section [5] by using the pseudo-differential calculus to 
establish Theorem 11.61 We have postponed the proof of Theorem 11.61 until this 
point as much of the needed notation will be established in Section [4] We will also 
complete the proof of Theorem 11.81 

We have chosen to use special case calculations, the functorial method, and the 
pseudo-differential calculus as our purpose in this paper is at least in part expository 
and we wish to illustrate the interplay amongst these methods. In a subsequent 
paper, we shall perform a similar analysis for other elliptic boundary conditions 
(Robin, transfer, transmittal, etc.); it will be necessary to restrict to Re(a) < 1 to 
ensure convergence and regularization will not be required in that analysis. 

2. Computations in M 2 

This section is devoted to the proof of Theorem 11.51 and we shall adopt the 
notation of that Theorem throughout. As wc shall be dealing with different weights, 
we drop the notation X and return to ordinary integrals in this section to perform 
a special case calculation in flat space. One has the following estimate of R. Lang 
[TO] and of H. R. Lerche and D. Sicgmund ^lTj that adjusts the formula of Equation 
()l.d|) for the heat kernel on a half space to take into account the curvature of the 
boundary of M (for related results see also [3J IT2"]): 

Lemma 2.1. Adopt the notation of Theorem \1.5i Let x = (y,r) g C £ . As t J. 0, 



4?rf 



PM (x, x; t) = -i- j 1 - e"'" 2 /' - L^^r 2 *" 1 / 2 H e^'dn } + 0(1) . 



rt-!/2 



Proof of Theorem \1.5\ (1). Parametrize the boundary of M by arclength. There 
is no higher order correction in R 2 and the 0(r 2 ) term in Equation (|l.fj) vanishes. 
Thus on the collar C e , we have 

dx = (1 - L aa (y)r)drdy . 

Following Equation (|l.d[) . wc let p H (x,x;t) = (47rt) _1 (l — e -1 " 2 /*) be the Dirichlet 
heat kernel in the halfspace r > on the diagonal. We take < a < 1 and express: 

Tr^Fe-*^') 

Fo{y)r~ a x{r){l - L aa (y)r)p M ((y, r), (y, r); t)drdy 



dM JO 

oo 

Fo(y)r " a xM(! ~ L aa (y)r)p M {{y, r), (y, r);t)drdy 

dM Jo 

D. + D2 + D 3 +D 4 + D 5 



where, motivated by Lemma |2.1| wc have: 

D 1 := J- / F(x)dx, 
4vrf J M 

1 p POO 

D 2 := --r- / / F (y)r- a e- r2 ^drdy, 

47rf Jom Jo 

A3 := ^ j j F (y)r- a (l - x (r))(l - L aa (y)r) e - r2 ^drdy, 



dM JO 
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i p poo 

L>4 := — / / F Q (y)L aa (y)r 1 - a e- r2 / t drdy, 

47rt JdM JO 

D 5 ~ [ f F (y)r- a x(r)(l-L aa (y)r){(p M -p H )((y,r),(y,r);t)}drdy. 

JdM JO 

A straightforward computation yields: 



We use Lemma [2TTI to compute D5. Since Fo(y)r a xi r ) is integrable on C e , we have 
the O(l) in Lemma [27X1 remains O(l) as t [ 0. Hence for < a < 1, 



1 />oo />oo 

D b = -— / / Fo^Laaiyy-ot-^e-^dvdrdy 

F (j/) J L aa ( 2 ;)L hb ( 2 ;)r 3 ~ Q <- 1 / 2 e - I ' 2 ^rfrd ? / + 0(1) 

1/2 

-r fc^) t- a ' 2 J^F {y)L aa (y)dy + O{y-^l 2 ) + O{l) 



1 

4?rf 



4tt 4(a-3) 
1 2 - a /2 - a 



r (^) ra ' 2 j dM F o(y) L ^y) d y + W 



4tt 4(a - 3) 

We collect terms to complete the proof of Theorem II .51 (1). □ 

Remark 2.2. We have chosen to study the region < a < 1. The reason for this 
is that F has to be integrable in order to control the 0(1) remainder in Lemma |2~T1 
If one wishes to obtain just the leading asymptotic behaviour of Tr L 2(Fe~ tAM ), 
then probabilistic estimates for 

(2.a) RAi(x;t) := {p M ~ PH){{y,r),{y,r);t) 

along the lines of pQ, and analogous to [2], could be used to show that for 1 < a < 3, 
(2.b) f f X (r)r' a R(x;t)drdy = 0(t- (m - 2+a ^ 2 ), i|0. 

JdM JO 



Proof of Theorem \1.5\ (2). Recall from Equation (jl.dj) the formula for the heat 
kernel pn on the diagonal for the half space. We decompose 

Tr L 2(Fe- tAM ) = Ex + E 2 + E 3 

where we have: 

E 3 ~ [ f F Q (y)r' 1 x(r){l~L aa (y)r){p M ^ PH ){{y,r),(y,r);t))drdy, 

JdM JO 

-1 /> />oo 

E 2 := --r- / / F (y)L aa (y) X (r)(l - e^^drdy, 
47rf J 9M Jo 

Ei:=— / ^ (2/)r- 1 x(r)(l - e^^drdy. 

4?rf JdM Jo 

We apply Equations (|2.a| and (|2.bj) with m = 2 and with a = 1 to see that 
E3 = Oy 1 / 2 ). Furthermore, E 2 = 0(t~ 1//2 ). The leading term is provided by E\. 
We integrate by parts to see that: 

Ei = ~ J gM Fo(y) j™ ln(r) | X '(r)(l - '*) + xWye" 1 " 2 /* J drdy 
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+0( e - £ o/(2*)) 

1 1 f nw. J, (4\ , „, n f e x(r) 



4?rf 2 



F (2/) In (#) + C + 2 / *^dr + O( e -o/( 2t )). 
This completes the proof of Theorem 11.51 □ 

3. The functorial method 

We adopt the notation of Theorem 11.61 throughout this section. We begin our 
study with the following: 

Lemma 3.1. There exist constants e v t a so that 

(1) a bd a (F,V) = (47r)-™/ 2 Z brf {Tr( £ Q F Id)}. 

(2) a\ d a {F,V) = (47r)- m / 2 X bd {Tr(e? iQ F 1 Id+£l iQ FoL aa Id)}. 

(3) a\ d a {F,V) = (4 7 r)-" 1 / 2 X bd {Tr(4 Q F 2 Id+£i iQ F 1 L QQ Id 

F E)}. 

Proof. We apply dimensional analysis - we shall suppose that a (£ Z for the moment. 
Let c > define a rescaling g c := c 2 g. We then have 

dx c = c m dx, dy c = c m ~ x dy, V c = c~ 2 V, 

r c :=cr, d Tc = c~ 1 d r , F l ^ c = c a ^ t F ll 

7" fWT fbd n m—l r rbd 

xRcg.c — c -^-Rcg7 -L-c — k x . 

Let a„. c := a n (x,T> c ), a e ^ c := a^ a (y,V c ), and ai <a ,i,c '■= a t,aA{y,'D c ) denote the 
local heat trace invariants defined by T> c on M and on dM, respectively. It is 
immediate that 

(3.a) Tr L 2(Fe- {tc ~ 2)v ) = Tr L 2 (Fe~ tv °) . 

We expand both sides of Equation (|3.a|) in an asymptotic expansion: 



t-^ 2 c m J2(c- 2 t) n X Kes {Fa n } 

n=0 

oo I 

+ {c- 2 t)-^- 1 ^ 2 Y J C a -^ l ^/ 2 Y J I bd {F l a l ^ i \ 



1/2 E 

n=0 



t n c m XR eg { Fa n t c } 



+t -(m-l)/2 t^' 2 J2 ^-^-^{FiOt^o}. 

e=o i=o 

Since a (£ Z, the interior and the boundary terms decouple. We equate terms in 
the asymptotic expansions to see that 

— 2n j i—l 

a n ,c = c a„ and a<i t a,i,c = c ai j0lt i . 

Examining relations of this kind is straightforward - they mean that the local 
formula a n (x, T>) is homogeneous of weighted degree 2n in the jets of the derivatives 
of the symbol of T> and that the local formula at ai i(y, T>) is homogeneous of weighted 
degree I — i in the jets of the derivatives of the symbol of T>. One may use Weyl's 
theory of invariants to express a spanning set for the invariants which arise in 
this context and complete the proof of Lemma 13.11 for a $ Z. We use analytic 
continuation to establish Lemma l3~H when a = 0, — 1, —2, ... as well. We refer to [7] 
for further details concerning this sort of dimensional analysis. 
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If a = 1, then the argument is rather different. Let e c be the width of the collar 
C with respect to the rescaled metric. The regularizing term in Equation ( |l.g[ ) docs 
not simply rescale. Rather we have: 

ln(e c )X b c d {(Fa n , c ) Q ,c} = hx{^)c a c- 2n c m - x l bd {{Fa n ) } 
= C " l - 2 "{ln(c) + He)}J bd {(Fa n ) } . 

This yields the modified relation: 

iRcsAFdnA = c m - 2n l Rcg {Fa n } + hx(c)c m - 2n l bd {(Fa n ) } . 

A similar argument for a = 2 shows that: 

ZReg,c{i^«,c} - c m - 2n I Rcs {Fa n } + lnfcjc™- 2 "^!^)! - (Fa n ) L aa } . 

When we use Equation (|3.a[) to equate coefficients in the asymptotic series, we have 

a bd a (F,V c ) = c e a bd a (F,V), 

which completes the proof of Lemma 13.11 in these exceptional cases. We compare 
the terms involving ln(c) to obtain additional relations. The left hand side in the 
following equation arises from ln(c _2 i) and the right hand side arises from T r eg,c 
when we apply Equation (|3.a[) : if k = 2n, then 

-21n(c)a&(F,X>) = Hc)l M {(Fa n ) Q } 
-2]n(c)a b k d 2 {F,V) = ]n(c)X bd {(Fa n )i - (Fa n ) L aa } . 

There are no corresponding terms if k is odd and thus a,k,i =0 and a,k,2 = if fc is 
odd. This establishes Theorem ll.6l (4). An alternate proof is given in Section 5. □ 

We now use the functorial method to establish the following result: 
Lemma 3.2. 

(!) 4,a = \ e l, a ande% a = 4 >a . 

(2) The constants of Lemma \3.1\ are dimension free. 

(3) £l, Q = £ 0,a-l; £ 2,a = £ 0,a-2> an( ^ £ 2,a = £ l,a-l- 

Proof. Suppose that M = Mi x A/ 2 , that gm = <?Mi + 9m 2 , that T>m = T^Mi + T^m 2 , 
and that Fm = F1F2 where Fj are defined on Mj. We suppose that Mi is a closed 
manifold and thus dM = Mi x dM%. We then have: 

g—tT>M _ g— tT>M 1 g — tT>M 2 

Tr L2 (F M e- WM ) = Tr^e"*^) • Tr^e" 41 ^) . 
Equating asymptotic series yields 

a bd a (F M ,V M )= a k (F u V Ml )a bd a (F 2 ,V M2 ) 

2k+j=l 

and hence a corresponding decomposition of the local formulas: 

(3.b) a b e d a l (y,V M ) = a k( x u T> Mi)a bd a , l {y2,T>M 2 ) for y = (x 1 ,y 2 ) . 

2k+j=t 

Assertion (1) now follows from Theorem 11.31 (2). from Lemma |3~T1 and from Equa- 
tion (|3.b|) ; the multiplicative constants (47r) _m / 2 play no role. If we take Mi = S 1 
and 2? Mi = — d 2 , then the structures are flat. Thus oq(x > 'Dm 1 ) = l/y/^n and 
afc(a;,"Djvfi) = for k > 1. Thus Equation (|3.b|) yields in this special case the 
following identity from which Assertion (2) follows after taking into account the 
multiplicative constants (47r)~ m / 2 : 

= ^a bd a /y 2 ,VM 2 ) ■ 
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Wc prove Assertion (3) by index shifting. Let x( r ) be a smooth function so that 
X = near r = e. Let F(y,r) = Fo(y)x(r)r~ ai for Rc(ai) < 3. Wc apply Theorem 
1.61 with a — ai and with a = a.\ — 1 to see that: 



4%,M V ) = 4-i, a -i,j-i(y,v) for j > i . 

Assertion (3) now follows. □ 

Proof of Lemma \1.9[ Assertion (1) of Lemma 11.91 follows from Lemma 13.11 and 
Lemma I3~2l by a suitable relabeling of the coefficients. Assertions (2) and (3) follow 
from Theorem 11.51 □ 

Proof of Theorem \l.£\ We now derive Theorem II .81 from Theorem 1 1 . 71 using Asser- 
tion (3) of Theorem 1.6. Certain of the coefficients are regular so the computation 
is elementary; we do not need to drop the pole. We simply set a = 1 to compute 
aW and a — 2 to compute a,Q d 2 and a| d 2 : 

a$ 2 (F,V) = K 2 (47r)^ m / 2 X M {Tr(- J F 1 Id)} 

= 0F(47r)- m / 2 X M {Tr(F o Id)}, 
a\ d x {F,V) = Ko (47r)- m / 2 X M {Tr(- J F 1 Id+|F L Qa Id)} 

= ^(4 7 r)—/ 2 2: M {Tr(-4F 1 Id+3F L aa Id)}, 
af a {F,V) = Ko (47r)-™/ 2 X bd {Tr(- J F 2 Id + 2 ^ T F 1 L aQ Id 

+ ^FoR amma Id — 3^ t _ 6 \ FoLaaLbb Id + 4(^_6) FoLabL a b Id 
~3(l-a) F oRijji l d-j^F E)}\ a=2 

= 0F(47r)-™/ 2 X M {Tr(-iF 2 Id +f F x L aa Id +±F R amma Id 
— ■§lF L aa Lbb Id +-g 2 -FoL a bL a b Id +^F i?y_,i Id +F n E)}. 

We compute the remaining coefficients as follows. Let c(a,i d a ,A) be the coefficient 
of the monomial I m {Ty(A)} in (^) m / 2 a\ d a (F,V). It follows by Theorem O that 

a^(F,P) = § • (4^)-" l / 2 X M {Tr(F Id)} . 

We may expand 

q-4 _ 2(q-3) + (2-q) _ . p ;4-a\ 1 

^a- 1 2(q-3) — 2(a-3) ~~ ^ L \ 2 ) 2(a-3) 1 

It now follows that 

c ( a i d 2i FoL aa Id) = — — 

a\%{F,V) = (47r)-™/ 2 X M {Tr(f Fild-(f + ±)F L aa Id)} . 
Many of the terms in a\ \ (F, T>) are in fact regular at a = 1 . We have: 
c(a b 2 d l7 F 2 ld) = -I c^i^ Id) = i, 

c ( a 24: Ramma Id) = ^2 ' c ( a 2Ji ^aa^bb Id) = — 49 ) 
c i a 2 d l> Lab L a b Id) = Jq . 

The terms involving -R^ji and -E can be written in the form —K a Fo(^Rijji Id+E). 
Thus we may use the regularization of —Fq Id in o,q \ which was already computed 
to see: 

agf^F.D) = (47r)-™/ 2 X bd {Tr(-iF 2 Id+iF 1 L aQ Id+^i? amma Id 

— ^L aa Lbb ld+-^L a bL a bId+j^Rijji Id+^E)} . 
This completes the derivation of Theorem 11.81 from Theorem 11.71 □ 
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4. The pseudo-differential calculus 

We adopt the following notational conventions. Let a = (pti, a m ) be a multi- 
index. We set: 

\a\ = a 1 + ... + a m , a\ 

t<S — a-x r r a m JO 

JU " KL ^ /\ . . . /\ iXJ jy^ J X 

We apologize in advance for the slight notational confusion involved with using a 
to control the growth of F and also to using a as a multi-index. 

We begin our discussion by reviewing the standard pseudo-differential com- 
putation of the resolvent on a closed manifold without boundary and refer to 
[3 [3 [H [13] for further details. Let V be an operator of Laplace type. We 
want to construct the resolvent of D — A for large A where we use Equation (|l.b[) to 
express: 

\&\<2 

For the symbol a(T>)(x, £) of D this means 

*(D)(x,0 = ]T a s (x)f; 

\S\<2 

note that for the scalar Laplacian Am the th term vanishes so ag = in this 
setting. In the evaluation of the heat equation asymptotics homogeneity properties 
of symbols are relevant and it turns out that collecting terms according to 

a 2 (x,£,\) = -A+ Y a s(x)C 

1*1=2 

a,j(x,€,\) = ^2a s (x)C, j = 0,l, 
|S|=j 

is fruitful. As a result, the symbol a(T> — \)(x,£, A) can be written as 

2 

a(V - \)(x,Z,X) =Y,ai(x,Z,*)- 

3=0 

For the symbol of the resolvent of T> — A we make the Ansatz 

oo 

(4.a) a((V - A)- 1 )^, £, A) ~ £ q-2-i(x, £, A). 

Z=0 

In view of the formula for the symbol of a product we see that q~ 2 -i is determined 
algebraically by 

(4.b) 1 = a 2 {x,t,X)q- 2 (x,t,X), 

(4.c) = V 4r[df%(x,e,A)] [U|g_ 2 _,(x i e,A)] for fc > 1. 
* — ' a! 

fe=2+i+|5|-j 

These symbols will play a crucial role in the proof of Theorem 11.61 in Section [5] 

We now specialize to the case where T> = Am is the scalar Laplacian. For our 
present considerations we need q- 2 , 9-3, q~A- For simplicity we skip the arguments 
in the following summary of results - we use in an essential fashion the fact that 
Am is scalar. We also use the convention that repeated indices are summed over. 



= a x \ x ... x a m l, 

= (jL-\ ai x X (-§J\ 

\ dxi J "' \ dx m J 
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Wc find (all greek indices will range over {1, 2, m}): 
-i 



9-2 

9-3 



"2 ' 
-1 



9-4 = ~a 2 1 



ai9_ 2 + (D%a 2 )(V-lD„q_ 2 ) 
a q-2 + aiq- 3 + (D^ai){\f^lDlq- 2 ) 
+(D%a 2 )(V=lDlq_ 3 ) - ±(D^a 2 )(D^q_ 2 ) 

For later use it will be advantageous to express the results in terms of q™ 2 . We then 
have 



9-3 
9-4 



-«l9-2 + C-3,39-2! 

-a q 2 _ 2 + C_4,39^2 + c -4,49-2 + C-4,59-2> 



where 

C-3,3 
C-4,3 
C-4,4 
C-4,5 



-V^l(^a 2 )(d», 

a 2 -V^T(d^ ai )(d x v a 2 ) - V=l(d%a 2 )(d* ai ) - ±(d^a 2 )(^a 2 ), 
-3a lC _ 3 ,3 + V^T(^a 2 )(^c_ 3 ,3) + (c^a 2 )(d> 2 )(d*a 2 ), 



= 3c^ 



3,3- 



The relevant operator for our considerations is 

A M - A = -g^d^d v + - A 

where wc have changed notation slightly from that used previously. For the symbols 
this gives 

a 2 (x,£,\) = g^^ u -\=\£\ 2 -\ 7 

ai(a;,£,A) = V^lb^^, and a Q (x, £, A) = . 

To state results for g_ 2 , q_ 3 and q_4 for this operator D we will as usual raise and 
lower indices using the inverse metric and the metric. Furthermore, "," denotes 
partial differentiation. One computes easily that: 

q- 2 (x, £, A) = npr=x 3 



q- 3 (x,£,X) = - m i_ x) i ' 
9-4(0;, C, A) = 



(ieP=A) 



+4g%g?Sg v %M>yZs + ^g^g^Up^} 



+7 



{-12g^g^g^g^U^s& } • 



'(l«P-A) 5 

If the manifold has a boundary the expansion (|4.aj) has to be augmented by a 
boundary correction. To formulate the conditions to be satisfied by the boundary 
correction we expand about r = 0. We adopt the notation established in Section 
11.71 and expand 

ds 2 M = 9tq(Vi ^dy" 7 ° dy e + dr 2 on C e . 
The coordinate y locally parameterizes the boundary, and r is the geodesic distance 
to the boundary, so x = (y, r). A tilde above any quantity will indicate that it is to 
be evaluated at the boundary, that is at r — 0. Furthermore, we use £ = (u>, t). 
Wc find 



A 



M 



k=0 \3\<2 



r=0 



V,r 
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with the notation 



i=l 



Introducing 

dj{y,r,uj,D r ,X) = 



for J = 0,1, 



E\ S \=2^(y^)[Ut~i 1 < i )D?--X for 3 = 2 
wc define the partial symbol 



As it turns out, the symbols 



fc! ^ dr> 



r=0 



2 00 1 d k 
aP>(y,r,u,D r ,\)=22 ^ —r k —^ai(y,r,uj,D r ,X) 



1=0 k=o 
l-k=j 



r=Q 



have suitable homogeneity properties and using these symbols we write 

2 

a'(A M -X)= J2 a {3) {y,r,u,D r ,X). 



j=-oo 



We write the symbol of the resolvent as 



(4.d) 



= ^2q-2-j(y,r,u,T, A) - e _ ^ TTr ^/i-2-j(2/,r,w,r, A), 

3=0 3=0 



where the second term is the boundary correction. The factor e ^ ^ TT appears be- 
cause the operator constructed from these terms is the Op'(h) in [13], and Op'(h) = 
Op{he-^ =lTr ). This shows 



a'(A M - A) o h-2-j (y, r,w,r, A) = 0. 

3=0 

Here o denotes the symbol product on H£ m_1 . Analogously to Equations (|4.b[) and 
(|4.c[) this equation leads to the differential equations 



= a (2) (y,r,w,D r ,A)/i_ 2 (y,r,a;,r, A), 
= a' 2 ) (y, r, u, D r , X)h-2-j (y, r, u, r, A) 

+ E 4 fe (fe W,w,£> r ,A)l f(V=lD w )a/i-a-i(l/,» , > w > T > A) 
' — ' a! L J L 

j=2+2+|S|-fc 
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For the present considerations we need h- 2 —j for j = 0,1,2, and we have more 
explicitly (repeated letters a, b, c, ... run over tangential coordinates {1,2, ...,m — 1}) 

= a (2 ' , (y,r,u,D r ,X)h- 2 (y,r,uj,T,X), 

= a (2) (y, r, u, D r , X)h_ 3 (y, r, u, r, A) + a (1) (y, r, u, D r , X)h_ 2 (y, r, u, r, A) 
+ \DlaPXy,r,u,D r ,\)\ [(V=TD v ) 6 ft_ 2 (i, J r > w J T > A)] , 

= a (2) (y, r, w, L> r , X)h- 4 {y, r, w, r, A) + a (0) (y, r, w, £> r , X)h- 2 (y, r, ui, r, A) 
+ [Djy i >( tf ,r,Ci/,A.,A)] [(\/=T£' I ') 6 /i_ 2 (2/,r,w,r,A) 

+| [lfca&>{y,r,u,D r ,\j\ [(^lD v ) bc h. 2 (y,r,uj,T,X) 

+a (1) (y, r, u), D r , X)h_ 3 (y, r,u),r, A) 

+ [D b y 2 \y,r,u,D r ,X)\ [{y/=lD v ) b h- 3 (y,r,u,T,X) . 

The relevant equations for a^- 1 (y, r, w, Z) r , A), i = 0, 1, 2 are 

a {2) (y,r,uj,D r ,X) = a 2 (y,r,u>,D r ,X)\ r=0 
= g ab cj a cj b + D 2 r - A, 

r(d r a 2 (y, r, w, D r , A))| r=0 + ai(y, r, u, D r , A)| r=0 
~Vb a oj a + \f^lb r D r , 



a^(y,r,u,D r ,X) 



rg ab uj a uj b + 

1 ^.2 / q2 



r (d r a 2 (y, r, u, D r , A))| r=0 + r{d r a x {y, r, w, D r , A))| r=0 
+ao(y,r,tjj,D r ,X)\ r=0 

±r 2 g ab r oj a uj b + ry^lb a r uj a + rV^lb r r D r + c. 



The differential equations have to be augmented by a growth condition 

(4.e) h- 2 -j(y, r, u>, r, A) — > as r — > oo, 

and an initial condition corresponding to the Dirichlet boundary condition 

(4i) h- 2 -j(y,r,(j,T, A)| r=0 = q- 2 -j(y,r,uj,T, A)| r=0 . 

Once the symbols /i-2-j have been determined, their contribution to the asymp- 
totics of the trace of the heat kernel follows from multiple integration. As before, 
we suppose r a F £ C°°(C E ). The contribution reads 



1=0 



Vi(y^ F , A M)dy 



DM 



with 



(4-g) 



^ (y ' F ' AM)= (^.?,L-^/ 



ds 



dre^ 



j+k=i • 
dre~ 



h-2-j {y, r , w, r, - V^T s)f fc - Q F fc (y) 



where 7 is anticlockwise enclosing the poles of h— 2 —j in the lower half-plane. The 
integral with respect to s is the contour integral transforming the resolvent to the 
heat kernel; see Section 5. Note that from ()4.dj) the contribution to the heat kernel 
is minus the above. 

As will become clear in the following, with A = 
of the type 



-Is, we need integrals 
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The r integration can be done using 



dre 



(\/^T) fc (-l) /+fe 7r / 1 d 



(r 2 + A 2 )' 



(I -I) 



2A dA 



1-1 



[A 



fe-l e -fAn 



So 



nkljn _ (v^T) fc (-l)"- 



ab... 



(' 00 

/ (is 



-1 s^l—a uj a u) b ... — rA /_1 d 



, fe-lg-fAl 



Performing the A-diffcrcntiation, different f-dependent functions would occur. It 
is therefore desirable to first perform the f-integration before performing the A- 
derivatives explicitly. This is achieved by noting that (z = A has to be put after 
the A differentiation has been performed) 

Hjn (x^T) fc (-l)" +fc+1 7 



T. 



(n-1)! 
1 d 



dw I dse s 



— s UJ a UJb- 



A" 



2 A dA 
(■ N /=l) fc (-l) n+fc+1 7r 



dff l - a e- f{K+z) 

z=A 

T{1 + 1 - a) / rfcj / dse^ 



AJ 



1 d 
2AdA 



A 



fe-i 



(A + z)^ 1 -" 

We can proceed in general by introducing numerical multipliers c n ki according to 

]_±\ n ^ A*- 1 
2A dA 



(A + z) 



l+l-a 



Cnki 



z=A 



1 

J^l+2n — k—a ' 



The s-integration is then performed using 

ds- 



2tt 



-M 2 



(l^p + V^Ts)/ 3 r(/3) 

The final ^-integrations follow from 



C(y) = 

by observing that 

_ ~ab 

dw Lo ai uj a . 2 ...uia r e 9 UJaUlb 
In particular 

dw e.-W 2 = 



1 V d 



8 



1 J dy ai dy a - 



C(y) 



7T 



'<)■ 



2 n 2 V99ab, 



du> Lu a LOi,e 

dw uj a u b uj c uj d e~ lujl = \fg{gab9cd + 9ac~9bd + 9ad9bc) 

Introducing the numerical multipliers d^iju according to 

2( v / ^T) fc (-l)" +fe+1 7r 2 r(/ + 1 - a)c nkl 



( n - l)!r ( j+ '- fe - a + n) 
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we obtain the compact-looking answers 

T ab. U . = d kl]n / duj CJ a CJ fc ...e" H 



where the last ^-integration is performed with the above results. 

Note that the numerical multipliers dkijn are easily determined using an algebraic 
computer program. Therefore, all appearing integrals can be very easily obtained. 

Let us apply this formalism explicitly to the leading orders, and we start with 
h—2(y,r,uj,T,X). The relevant differential equation reads 

(d 2 r -A 2 )h_ 2 (y,r,Lu,T,\) = Q, 

which has the general solution 

h- a (y,r,w,T, A) = Ae~ rA + Be rA . 

The asymptotic condition (|4.c|) on the symbol as r — > oo imposes B = 0. The initial 
condition h- 2 \ r =o = <?-2|r=o gives 

A=- 1 



, 2 + A 2 ' 

Putting the information together we have obtained 

h-2(y,r,uj,T,X) = r2 l +k2 e- rK . 
Performing the relevant integrals, with the notation 

oo 



s / dfe^ s \ - J (ire- 



produces 



dIh- 2 {y,r,w,T,-y/-ls) = doooif 2 V 9 = y (l - -) n 2 

Taking into account the prefactor in Q4.g[ ), this confirms the value of R a in Lemma 
Ol 

In the next order we obtain 

(d 2 - A 2 )/i_ 3 (y, r, w, r, A) = (E + U 1 )e~ rX + (F + U 2 )re~ rA , 

where 

b r A = gfwgUb 

r 2 + A 2 ' ~ r 2 + A 2 ' 

~ac. ,6. , , , / T 7,ac, ,b. 



C/1(W) - r 2 + A 2 + (r 2 +A 2 ) 2 ' ^ ~ (r 2 +A 2 )A ' 

Note, for later arguments, that U\{lj) and U 2 {ui) are odd functions in w. Fur- 
thermore, for the scalar Laplacian at hand b a = g hc T b c a \ thus they contain only 
tangential derivatives of the metric. 

Using for example the annihilator method, we write down the general form of 
the solution to this differential equation as 

h- 3 (y, r, u, t, A) = Cl e~ rA + c 2 re~ rA + c 3 r 2 e - rA + c 4 e rA . 



From the asymptotic condition (|4.ej) we conclude C4 = 0. From the initial condition 
given in Equation (|4.fj) we obtain 



Cl 



( T 2 +A 2)2 ( T 2 +A 2)2 ( T 2 +A 2)3 ( T 2 +A 2)3 
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From the differential equation we derive 

Collecting the available information, we see 

h- 3 (y, r, u, r, A) = Der rK + Br e - rA + Cr 2 e - rA + 0(w), 



with 



f-Lb r T 2y/-l g ab TLO a U b 



D = 
B = 
C = 



■ 2 +A 2)2 (r 2 +A 2 ) 3 

g ab io a uj b b r 



4A 2 (r 2 +A 2 ) 2(r 2 +A 2 )' 
g ab u; a uj b 



4A(r 2 + A 2 ) ' 



and where 0(lo) is an odd function in cu. Furthermore, 0{uj) contains only tangen- 
tial derivatives of the metric. We next perform the multiple integrals; note, odd 
functions in uj do not contribute. We obtain 



dlh- 3 (y, r, u>, r, s) 

n 2 \Zll9 a r9ab {— "^-^^1002 + 1^0101 — V - 1 ^1003 — ^0121 _ 1^0211 j 

7r(a - 4) (2 - a" 



4(3- a) 



„m/2 /^xab~ 
n \J99, r 9ab- 



This confirms the value of in Lemma 11.91 after taking into account the prefactor 
in Q and the fact that gfg ab = -g ab g a b,r = 2g ab L ab . 

Up to this point the calculation can be considered a warm up for the next order. 
We would like to determine the universal coefficients of the geometric invariants 
Ramma, L aa L bb and L ab L ab . In terms of the metric the last two are determined by 

T — — 1 n 
-L^ab — 29db,r • 

Using the Christoffel symbols 

r jfc* = (9lj,k + 9kt,j - 9jk,l) , 

and taking into account that with our sign convention the scalar curvature is given 
by the contraction g^Rijk 1 , we may expand the Riemann curvature tensor in the 
form: 

p l I I I "P ^F 1 n P ^F 1 n 

*Hjk — ^ jk .1 *■ ik ,j ~r 1 in jk ^ jn ik 

The normal projection of the Riemann curvature tensor reads 

£> Ln ac n - Ln ac n - l-n^c-ad- ~ 

^ararna 29,r9ac,r 2.9 9ac,rr $.9 ,r 9 : r 9caybd 

l_n a b~cd~ ~ l~OC~ 

— 4^9 9 yac,rybd,r 2" yac,rr- 

The above results suggest a strategy for the calculation. It suffices to consider the 
special case where the metric is independent of y. As a consequence, our answer 
will have the form 

(4.h) (4ir)- m / 2 {Ag ac g a c,rr + Bg ab 9 cd 9ac,r9bd,r + Cg ab 9 cd 9ab,r9cd,r} ■ 

This has to be compared with the terms in a b 2 a {F, Am) that possibly contribute 
to these geometric invariants. In detail one can show these terms are (mod terms 
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with tangential derivatives of the metric) 

hn a R + K^yRamma + K n Laa,Lbb + l^,L a bL a 



gK Q rf, -f- K a rL a mma T K a -L aa^bb T ri a Li a bl^ab — 

-ac- H„ 1 „.3 \ i -a6-cd- - / 1 „ , lj 1 ! 

y yac,rr \g h a J a / + " » yab,rycd,r ^24 K a T 4 K aJ 

~ab~cd~ ~ I 1 _i_l„3 , 1„5' 



So once we know A, £?, C, we can deduce 

(4.i) «a =-2(A-i« a ), «*=4(<7-£« a ), «i = 4(B + i« a -i/4) . 

In summary, when writing down the differential equation for h—^(y, r, w,t, A), we 
can neglect all terms that are odd in w as well as all terms that contain tangential 
derivatives of the metric. We obtain (up to irrelevant terms) 

{d 2 - A 2 )/ l _ 4 (y, r, u, r, A) = Ae~ rA + Br £ - rA + Cr 2 e- rA + Dr 3 e~ rA : 

where 



A = 

B = 

C = 

D = 



/^Tb r fc r Ar , ^V^^g^ Atoj.hh _ b r g a r ui a u b __b^_ 

(t 2 +A 2 ) 2 p^+iV 5 ) 3 4A 2 (r 2 +A 2 ) + 2(r 2 +A 2 ) ' 

b r r A \f~^Lg a ^b r Tu) a u] b 2^19^9^™^^^ g° r b r w a ui b b r b r A 

~ t 2 '+A 2 (t 2 +A 2 ) 2 (V 2 +A 2 ) a 4A(r 2 +A 2 ) _ 2(r 2 +A 2 ) ' 



2(r 2 +A 2 ) 4A 2 (r 2 +A 2 ) ^ 4(r 2 +A 2 ) 



9 T g r ^a^b^c^d 



4A(t 2 +A 2 ) 

So the solution has the form, taking into account the asymptotic behavior (|4.eD . 

h_ 4 (y,r,w,r,A) = ae" rA + /3re" rA + 7 r 2 e~ rA + (5r 3 e" rA + er 4 e - rA . 
From the initial condition a = (7-4(2/, r, u>, r, A)| r= o we obtain, up to irrelevant terms, 

+ (t *+a*)* {-Qb r g a r b T 2 uj a LJb + 4gf r r 2 uj a uJb + 2gfg^u) a u) b w c uj d } 

+ ( T i+ A iy-, {-l2gfgfuj a uJbUJ c ujdT 2 } . 

From the differential equation we obtain the conditions 

A = -2A/3 + 27, £ = -4A 7 + 6(5, 
C7 = -GAS + 12e, L> = -8eA . 



This determines the numerical multipliers fi, 7, (5 and e to be 

= 



3D_1C_1B_ 1 A 
'8F 4F 4 A 2 " 2 A ' 



— 3D_1C_ IB 
' 8F 4 A 2 " 4 X' 

" 4. A 7 6 A ' 

fc 8 A • 

For the Laplacian on the manifold M we have 

— — 2ff ffab.r, 

i>~ ~ab 1 ~ab~cd~ ~ 

5ab.9,r — 2-9 .9 <7a6,r Scd.n 



6 r = 



2- 1 

1 ~ac~bd~ ~ 1 ~ab~ 



gabg^r = z 9"~ V" g a b, r gcd, r - g"~ g a b, rr , 
ff^fr (5ah.9crf + g ac gbd + gadgbc) = g ah g cd g a b. x gcd,r + 2g ab g cd g ac ,rgbd,r- 
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Performing the integrations we obtain the contributions (modulo tt^" 1 l " 2 y/g) 



g ab g cd g ab r g cd r [— |d 2 003 + ^0003 — §^2004 + f^0004 
~ab ~cd ~ 



3d 



+ 9 a 9° 9ac, r gbd,r [^2003 — <^0003 + 4d 2 004 + d 
+ 9° 9ab,rr [~^2003 + \d 



0004 



6d 



2005J 



0003 



2d 2 



g ab g cd g abr g cdr 



-d u 



+ g ab g cd gac, r gbd, r 



128^0151 + ^4^0151 



32^0111 



23 



16 



2005 



122 



103 



^^0151 + ^^0151 + —4-^^x123 — S^0131 + 1^0111 



32 



123 — 8«0131 



+ g ab gab,rr [^^0131 ~ g^Olll] 1 



11 



. 16 

~ab~cd~ ~ 

g g g a b, r gcd,r 



158^0241 + (=1^0241 



— r^-^1213 + -55^0201 — (^^0221 



til' 



+ "^^^1212 + £^0221 



1 ~ab~cd~ ~ 
+ g g 9ac,r9bd,r 

1 7.ab ~ 



g4<^0241 + ^2<io241 



"4—^1213 — 1^0221 + 1^0201 



g ab gab,rr [j^^0221 — |^020l] 



<>I = g ab g Cd gab, r gcd,r [jf^ 0331 ~~ 3^^031l] 

+ g ab g cd gac, r gbd,r [ggdo33i — tk^ 0311 ] 

+ g ab gab,rr [^do31l] ! 
£/ = g ab g Cd gab,rg C d,r [j5gdo42l] + 9°" 9° 9ac,r9bd,r [^<^042l] • 

Adding up all terms and simplifying using the functional equation and the doubling 
formula for the T-function, the contribution to the heat kernel coefficient reads 



fA--\— m/2 J ~ab~cd~ ~ 

(47T) j_g g g a b,r9cd,- 



3a -16a-27 p I !-« ' 
384(a-6) \ 2 t 



txab-cd- % 5(9+4a-a-) r (l-a\ *%ab%, a+3 r n-aU 
+9 g gac, r gbd,r i 9 2( a _6) 1 \2^> + g gab,rr^-L ))■ 



This allows us to read off A, B, and C from Equations (|4.h[) and (|4.i[) and to 
conclude: 



,_3 



Ua-lK, k 4 c 



12 



7-8a+a" 
16(a-6) 



8 (a- 6) 



This completes the proof of Lemma 11.101 

5. The proof of Theorem 11.61 



We adopt the notation of Theorem 11.61 We may suppose that the weighting 
function F is supported in a boundary coordinate neighborhood U and that F has 
the form F = r~ a G with G € C^ mp (U). For appropriate functions q and q bd set 



Op(q)f(y,r) := 
Op'(q bd )f(y,r) 



-ly-u-\-y/ — lrr 



<l(y,r,u,T, X)f(uj,T)dujdT, 



-ly-uj bd 



9 (y,r,u,T,\)f(u,T)dwdT, 



du> 



(2 7 r) 1 - m dw, and dr := (2n)- 1 d T . 



Thus one has that Op'(q bd ) = Op(e v^rTgbdy j n \ oca \ coordinates (y, r), we follow 
the construction in [13) to define the standard parametrix for (T> — A) -1 : 



N 



QnW = £ {PP(?-2-n(A)) - <V(^ a _ n (A))} . 



n=0 
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The qj and q- are discussed in further detail in Section |3J we refer in particular to 
Equations P~aj) . (|4.b|) and (gx]). Note that the q bd here is the hj of Section gj 

We pass to a parametrix Hm for the operator e~ tT> by performing a contour 
integration: 

/M2?,t) := — = / e- tx Q N {\)d\ = V{Op(e_ 2 _„) - Op'(e bd 2 _ n )} ; 

in this expression, ej and e bd are the corresponding contour integrals of qj and (7* . 
For another construction of the heat parametrix see [8] . 

Take 4> and tp with compact support in U, with ip = 1 in a neighborhood of the 
support of 0. Then 

fljv - [e~ tv - Hjv] V 

has a kernel k(y,r,y' ,r';t) which, for sufficiently large N, is C 2 , 0(t J ) for large J 
(depending on N), and vanishes when r = 0, although not for r' = 0. To achieve 
vanishing in both variables, we consider 

4>H N (V,t/2)i; 2 H N (V*,t/2)*(f>~0e- t ' D / 2 ip 2 (e- t ' D y2)* < p 
= -R N i>H N (V*,t/2)*<t> - <f>H N (p, t/2)TpR% - R N R* N . 

Here we compute the adjoints with respect to the measure dydr. The remainder 
given by Equation (|5.a|) has a kernel which is C 2 , with derivatives 0(t J ) for large 
J, and vanishes when either r = or r' = 0. By the pseudo-local properties of Hn 
and e - * 23 , the same is true of 

4>H N {V, t/2)H N (V*,t/2)*(j> - 4>e- w 4> . 

Hence the kernel of this operator is 0(r 2 t J ) for large J and we can deduce an 
expansion for Tr L 2(Fe~ tT> ) from an expansion of Tr^^FHNH^cf)), taking <j> = 1 
on the support of F. The error will be holomorphic in a for Rc(a) < 3. Since 
Tr^2 (Fe~ tT> ) is also holomorphic for Rc(a) < 3, we may compute our expansion for 
Re(a) < 1 and continue analytically. 

From the expansion of qj in powers of (\uj\ 2 + t 2 — A) -1 (See Section H|) we get 
the estimates 

(5.b) |e_ 2 _ fc | < C t k ' 2 e- Ct{ ^ 2+T2) 

(5.c) \ef\ < C(\uj\ 2 +t 2 + i)-i e "l-l 2t /2 

for suitably chosen constants C and c. Set £ = (w,r). The qj and the q bd have an 
appropriate homogeneity property [13j . This homogeneity yields that 

(5.d) ej(y,r,s£,t/s 2 ) = s 2+J e J (y, r, £, t), and 

(5.e) ef{y, r/s, s£, t/s 2 ) = s 2 ^ef{y, r, £, t) . 

Furthermore, the kernel k\ of Op{e.j)Op(e* k )* and the kernel ki of Op(ej)Op' (e k d *)* 
on the diagonal are given, respectively, by: 



(5.f) h(y,r,y,r;t) = j e 3 (y, r, £, t/2)e k (y, r, £, t/2)d£ 

= {t/2) -(m +j+ k+v/2 j ej ( y , r} t:,l)e k (y,r,Z,l)dZ, 
(5.g) k 2 (y,r, y> r;t) = ( t / 2 )-(«+i+*+4)/a / ej (y, r, £, l)e^(y, r(2A) 1/2 , £, . 



There are similar formulas for the kernel of Op' \& d )Op{e* k )* and for the kernel of 
Op' {e bd )Op' {e bd *)* on the diagonal. 

We integrate F(y,r)ki(y,r,y,r;t). An expansion of Ftjek in powers r l ~ a gives 
terms i-( m +J+ fc + 4 )/ 2 times a meromorphic function with simple poles at a = 1, 2, .... 
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For other a, this gives the interior terms in Theorem 1 1.61 the terms with j + k odd 
vanish because of the parity of ej and in £. For the terms in Equation ( |5.g[ ), note 
that from Equations ()5.c[) and (|5.ep , e| decays exponentially as r — > oo so that we 
may integrate in r from to oo. An expansion of Fej in powers of r and a change 
of variable r/^/i — > r gives boundary terms of the form in Theorem II .61 for a ^ 1,2. 

This proves Theorem 1 1 . 61 for Re(a) < 1 and the rest follows by analytic continua- 
tion. Since for each t, the expansion is holomorphic in a for Rc(a) < 3, the residues 
arising from the interior integrals at a = 1,2 must be cancelled by residues from the 
boundary terms. And, since the various powers of t are linearly independent, the 
residue of each interior coefficient must be cancelled by the residue from the corre- 
sponding boundary coefficient. The residue from t n In cg {F a n ) = t n I-R CS (r~ a Ga n ) 
is 

/ -t n l bd {{Fa n )o} at a = l, 

[ ' 1 -f2 M {(Fo n ) 1 + (Pa B )oU at a = 2. 

The corresponding boundary term when a = 1 is i" _ ( a_1 )/ 2 a2^ a (F,T>). Let R be 
the residue of a^„ a (F,T)) at a = 1. Then 

t n-{ a -l)/2 a bd jp ty = _^ _ t n ln ^ R + ( const)t n + _ y 

a — 1 2 

Assertion (2) now follows. Since t n R must cancel the residue in Equation (|5.hp . wc 
get the ln(t) coefficient in Assertion (4) for a = 1. The proof of the rest of Assertion 
(4) and also of Assertion (3) follows in a similar fashion. □ 
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